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PREFACE 


The  work  described  in  this  report  was  conducted  between  April  of 
1984  and  July  1985  by  personnel  in  the  Radar  Division  of  the  Environ¬ 
mental  Research  Institute  of  Michigan  (ERIM).  The  work  was  supported 
by  Intera  Environmental  Consultants,  Inc.,  as  part  of  U.S.  Geological 
Survey  Contract  No.  14-08-0001-21748.  The  technical  monitors  for 
this  program  were  Messrs.  Joseph  Pearson  (Intera)  and  John  Jones 
(USGS). 

The  principal  investigator  for  this  program  was  Mr.  Eric  S. 
Kasischke.  Mr.  Richard  W.  Larson  and  Dr.  Philip  L.  Jackson  assisted 
in  the  design  and  deployment  of  the  calibration  and  geometric  refer¬ 
ence  reflectors.  Messrs.  James  Marks,  Jack  Losee  and  Robert  Williams 
processed  the  SAR  data. 
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USGS  X-,  C-  AND  L-BAND  SAR  DATA  COLLECTION  PROGRAM 


INTRODUCTION 

This  report  presents  a  summary  of  the  multi-channel  ( i . e . ,  multi¬ 
frequency  and  multipolarization)  synthetic  aperture  radar  (SAR)  data 
set  collected  for  the  United  States  Geological  Survey  (USGS)  during 
the  Side-Looking  Airborne  Radar  (SLAR)  Data  Acquisition  Program. 
This  report  specifically  discusses  data  collected  over  the  North 
Carolina  Digital  Project  Area  (NCDPA).  This  SAR  data  collection  was 
conducted  between  8  and  10  April  1984  over  four  forested  test  sites 
in  central  North  Carolina  (see  Figures  1  through  3). 

The  objective  of  this  data  collection  program  was  to  provide  USGS 
scientists  with  a  data  set  which  could  be  used  to  determine  whether 
or  not  geologically-stressed,  vegetated  areas  could  be  discriminated 
using  multi-channel  SAR  data.  In  order  to  achieve  this  goal,  the 
SAR  data  must  have  a  high  degree  of  radiometric  and  geometric  fidel¬ 
ity.  To  assist  in  obtaining  this  fidelity,  a  set  of  calibrated  cor¬ 
ner  reflectors  and  a  set  of  geometric  reference  reflectors  were  de¬ 
ployed  near  the  four  test  sites. 

During  the  data  collection  period,  three  SAR  missions  were  flown, 
with  22  passes  of  imagery  being  collected.  Two  missions  (16  passes) 
were  flown  with  the  SAR  collecting  X-  and  L-band  data,  and  one  mis¬ 
sion  (6  passes)  was  flown  with  the  SAR  collecting  X-  and  C-band  data. 

There  are  five  chapters  in  this  report,  including  this  introduc¬ 
tion.  Chapter  2  describes  the  ERIM/CCRS  CV-580  SAR  System  used  to 
collect  the  data  for  this  project.  Chapter  3  documents  the  data  col¬ 
lection  and  processing  performed,  including:  (1)  coverage  maps  and 
SAR  system  parameters;  (2)  a  discussion  of  the  calibrated  corner  re¬ 
flector  array;  (3)  a  discussion  of  the  geometric  reflectors;  (4)  a 
discussion  of  the  digitally-recorded  data;  and  (5)  an  engineering 
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Figure  2.  Locations  of  Test  Sites  A,  B  and  C  of  the  NCDPA 
(after  USGS  Raleigh  1 :250000  Quad  Sheet) 
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Figure  3.  Locations  of  Test  Sites  D  and  E  of  the  NCDPA  (after  USGS 
1 .250000  Greensboro  Quad  Sheet) 
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assessment  of  the  SAR  data.  Finally,  Chapter  4  summarizes  the  prin¬ 
cipal  results  of  this  program. 

There  are  four  appendices  to  this  report.  Appendix  A  discusses 
the  basic  manner  in  which  the  X-  and  L-band  SAR  imagery  collected 
during  this  experiment  may  be  calibrated.  Appendix  B  describes  the 
dihedral  reflectors  deployed  during  the  SAR  data  collection.  Appen¬ 
dix  C  describes  the  Hybrid  Image  Processing  Facility  used  to  digitize 
the  SAR  data  for  this  program.  Appendix  0  discusses  the  formats  of 
the  digitally-recorded  SAR  data  set. 
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CV-580  SAR  SYSTEM  DESCRIPTION 


From  1978  through  1984,  ERIM  operated  a  mu  1 ti frequency  and  multi¬ 
polarization  SAR  in  conjunction  with  the  Canada  Center  for  Remote 
Sensing  (CCRS).  This  SAR  was  maintained  through  a  leasing  agreement 
between  ERIM  and  CCRS  and  is  called  the  CV-580  SAR  System.  This  fa¬ 
cility  was  used  extensively  over  the  past  six  years  in  a  variety  of 
oceanographic  and  terrestrial  research  programs.  The  leasing  agree¬ 
ment  between  ERIM  expired  the  end  of  CY84.  At  this  time,  the  CV-580 
SAR  System  is  not  available  for  U.S.  Government  sponsored  research 
programs.  Efforts  are  underway  to  transfer  this  system  into  a  U.S.- 
owned  aircraft  (Shuchman,  et  al . ,  1984). 

The  CV-580  SAR  System  was  essentially  an  X-band  radar  to  which 
an  L-band  and  a  C-band  capability  was  added.  Because  of  this  config¬ 
uration,  the  X-band  channels  were  always  available,  and,  in  addition, 
one  could  operate  either  the  L-band  or  the  C-band  channels  simulta¬ 
neously  with  the  X-band  channel.  For  each  frequency  band,  two  ortho¬ 
gonal  polarizations  were  available.  The  radar  could  transmit  either 
horizontal  or  vertical  polarization  and  receive  both  the  parallel 
and  perpendicular-polarized  returns.  The  radar  antenna  could  be 
pointed  so  as  to  image  on  either  side  of  the  aircraft.  Table  1  sum¬ 
marizes  the  relevant  radar  parameters  for  this  system. 

Figure  4  presents  a  functional  diagram  of  the  CV-580  SAR  System. 
The  major  components  of  the  SAR  are  its  antennas,  the  transmitters 
and  receivers,  the  calibration  signal  generator,  the  optical  film 
recorders,  the  digital  tape  recorder  and  the  X-band  real-time 
processor. 

Figure  5  presents  a  more  detailed  diagram  of  the  reference  oscil¬ 
lators,  transmitters  and  receivers  of  the  SAR  system.  This  diagram 
shows  how  the  L-band  and  C-band  wavelengths  are  returned  to  X-band. 


TABLE  1 

CV-580  SAR  SYSTEM  PARAMETERS 
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Figure  5.  Schematic  Diagram  of  CV-580  SAR  System  Transmitter/Receiver 
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Figure  6  presents  a  diagram  of  the  calibration  signal  generator 
(CSG)  system,  which  is  used  as  a  monitoring  and  reference  source  for 
the  periodic  verification  of  the  SAR's  system  transfer  function. 
The  calibration  signal  generator  produces  synthetic  target  signals 
which  are  added  into  the  radar  receiver  (Walker  and  Larson,  1981; 
Larson,  et  al . ,  1981;  1982;  1985).  The  form  of  the  synthetic  target 
corresponds  to  a  given  range,  and  the  intensity  is  controlled  by  an 
accurate  calibrated  R.F.  attenuator.  At  specific  times  during  the 
imaging  flight,  the  calibration  signals  are  inserted  into  the  re¬ 
ceiver  at  the  antenna  terminals,  so  that  the  radar  receiver  detects 
them  as  a  radar  return.  Using  a  sequence  of  such  signals  wi^h  dif¬ 
ferent  intensities,  along  with  the  processed  signals  from  the  corner 
reflectors,  a  calibration  curve  can  be  produced  of  output  signal 
power  versus  radar  cross  section  (a  ).  When  combined  with  measure¬ 
ments  of  precision  corner  reflectors,  CSG  measurements  can  be  used 
to  achieve  absolute  calibration  of  the  SAR  (see  Walker  and  Larson, 
1981;  or  Larson,  et  al.,  1985).  The  synthetic  targets  from  the  CSG 
can  be  recorded  both  optically  and  digitally,  but  only  at  X-  and 
L-band. 

The  three  SARs  which  comprised  the  CV-580  SAR  System  are  similar 
in  that  they  all  use  the  synthetic  aperture  technique  to  produce 
imagery  with  fine  cross-range  resolution.  These  SARs  used  pulse  com¬ 
pression  to  achieve  fine  resolution  in  the  range  dimension.  Each 
recording  channel  was  adjusted  to  produce  imagery  of  a  selected  swath 
parallel  to  the  flight  direction.  The  width  of  the  imaged  swath  was 
determined  by  the  range  increment  of  signals  displayed  on  the  record¬ 
ing  system.  The  displacement  of  the  recorded  swath  from  the  flight 
line  was  adjustable  by  the  radar  operator.  The  CV-580  SAR  System  is 
extensively  described  by  Rawson,  et  al.  (1975). 

The  SAR  phase  histories  were  recorded  onboard  the  aircraft  both 
optically  on  70  mm  film  and  digitally  on  high  density  digital  tape 
(HOOT).  The  optical  recording  system  recorded  four  data  channels, 
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while  the  digital  system  was  restricted  to  recording  only  two  of  the 
four  data  channels  during  data  collection.  In  addition  to  the  opti¬ 
cal  and  digital  recorders,  a  real-time  digital  SAR  processor  onboard 
the  CV-580  generated  X-band  imagery.  The  real-time  imagery  was  not 
intended  to  be  of  high  quality,  but  provided  "quick-look"  data  neces¬ 
sary  for  system  performance  evaluation  and  planning  purposes  during 
data  collection  missions.  It  was  also  used  to  select  digital  data 
of  interest  for  subsequent  processing  and  analysis.  Table  2  presents 
the  various  data  recording  options  available  on  the  system,  figure 
7  schematical ly  illustrates  the  various  SAR  data  recording  and  pro¬ 
cessing  options. 
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This  chapter  contains  a  summary  of  all  multi-channel  SAR  data 
collection  missions  conducted  during  the  USGS  North  Carolina  program, 
as  well  as  a  description  of  all  SAR  data  processing  performed  at 
ERIM.  This  chapter  also  contains  a  description  of  the  calibration 
and  geometric  reflector  arrays  which  were  deployed,  as  well  as  a 
description  of  the  digital  SAR  image  products  which  were  generated. 

3.1  SAR  DATA 

There  were  five  test  sites  within  the  North  Carolina  Digital 
Project  Area  (NCDPA).  These  sites  are: 

A.  Cedar  Rock/Coolers  Rock  Mountains 

8.  Shepard/Caraway  Mountains 

C.  Pilot  Mountains 

D.  Bowlings  Mountain 

E.  Daniels  Mountain 

The  locations  of  these  sites  are  presented  in  Figures  1  through  3. 
Three  SAR  data  collection  missions  were  flown  over  these  sites  on  3 
through  10  April  1984  with  the  CV-580  SAR  System.  A  total  of  22 
passes  of  SAR  imagery  were  collected.  These  missions  are  summarized 
in  Table  3. 

During  the  NCPDA  data  collection,  the  CV-580  SAR  System  was  con¬ 
figured  in  its  two-frequency,  narrow  swath  mode.  In  this  mode,  four 
channels  of  SAR  data  were  recorded  optically  on  SAR  signal  film  and 
two  channels  recorded  digitally  on  high  density  digital  tape.  Each 
optically-recorded  channel  has  a  slant-range  swath  width  of  5.4  km 
and  each  digital  channel  has  a  slant-range  swath  width  of  6.1  km. 

In  the  two  frequency,  narrow  swath  mode,  either  horizontally  (H) 
or  vertically  (V)  polarized  microwave  energy  is  transmitted  in  two 
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frequencies  (X-  and  L-band  or  X-  and  C-band).  Both  like  (HH  or  VV) 
and  cross  (HV  or  VH)  backscattered  energy  from  each  frequency  is 
recorded  optically,  resulting  in  four  channels  of  data.  Two  of  the 
four  channels  are  recorded  digitally. 

X-  and  L-band  SAR  data  were  collected  on  8  and  9  April  1984  and 
X-  and  C-band  data  were  collected  on  10  April  1984.  The  frequency 
and  polarization  combinations  which  were  optically-  and  digitally- 
recorded  are  also  summarized  in  Table  3. 

In  the  remainder  of  this  section,  we  will  discuss  the  SAR  data 
collected  during  each  of  the  three  missions.  For  each  SAR  mission, 
a  table  summarizing  the  aircraft  ground  tracks  and  SAR  parameters 
are  presented.  Diagrams  which  present  t1";  ground  location  of  each 
pass  are  also  included,  as  are  representative  examples  of  SAR 
imagery. 
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3.1.1  USGS-1  (8  April  1984  -  13:06  to  16:24  EST) 

The  purpose  of  this  mission  was  to  collect  multipolarization 
X-  and  L-band  SAR  data  over  the  NCDPA.  Three  passes  of  SAR  imagery 
were  collected  over  a  test  area  to  verify  the  operational  status  of 
the  CV- 580  SAR  System.  Six  passes  of  SAR  imagery  were  collected  over 
Test  Sites  A,  and  B  and  C,  including  one  pass  over  the  calibrated 
corner  reflector  array. 

The  aircraft  and  SAR  system  parameters  for  the  USGS-1  mission 
are  summarized  in  Table  4.  In  Table  4,  all  latitudes  and  longitudes 
represent  the  positions  to  the  near  edge  of  the  swath.  The  locations 
of  the  ground  coverage  during  each  SAR  pass  of  the  USGS-1  mission 
are  presented  in  Figures  8  through  15.  Figure  16  presents  optically- 
processed  X-  and  L-band  SAR  imagery  of  Pilot  Mountain  collected  dur¬ 
ing  Pass  7.  The  SAR  data  collected  during  this  mission  was  all  of 
excellent  qual i ty. 


SUMMARY  OF  SAR  AND  AIRCRAFT  PARAMETERS  DURING  USGS- 
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3.1.2  USGS-2  (9  April  1984  -  11:27  to  13:45  EST) 

The  purpose  of  this  mission  was  to  collect  multipolarization 
X-  and  L-band  SAR  data  over  the  NCDPA.  Two  passes  of  SAR  imagery 
were  collected  over  a  test  area  to  verify  the  operational  status  of 
the  CV-580  SAR  system.  Five  passes  of  SAR  imagery  were  collected 
over  Test  Sites  A,  8,  C,  0,  and  E  including  one  pass  over  the  cali¬ 
brated  corner  reflector  array. 

The  aircraft  and  SAR  system  parameters  for  the  (JSGS-2  mission 
are  summarized  in  Table  5.  The  locations  of  the  ground  coverage  dur¬ 
ing  each  SAR  pass  of  the  USGS-2  mission  are  presented  in  Figures  17 
through  21.  The  SAR  data  collected  during  this  mission  was  all  of 
high  quality. 
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3.1.3  USGS-3  (10  April  1984  -  15:59  to  17:58  EST) 

The  purpose  of  this  mission  was  to  collect  mul ti polari zation 
X-  and  C-band  SAR  imagery  over  the  NCDPA.  One  pass  of  SAR  imagery 
was  collected  over  a  test  area  to  verify  the  operational  status  of 
the  CV-580  SAR  System.  Five  passes  of  SAR  imagery  were  collected 
over  Test  Sites  A,  B,  C,  0,  and  E  including  one  pass  over  the  cali¬ 
brated  corner  reflector  array. 

The  aircraft  and  SAR  system  parameters  for  the  USGS-3  mission 
are  summarized  in  Table  6.  The  locations  of  the  ground  coverage  of 
each  pass  of  the  USGS-3  mission  are  presented  in  Figures  22  through 
25.  Figure  26  presents  optically-processed  SAR  imagery  of  Pilot 
Mountain  collected  during  Pass  4.  The  SAR  data  collected  during  this 
mission  was  all  of  excellent  quality. 


I  ABIE  6  SUMMARY  OF  SAP  AND  AIRCRAFT  PARAME1ERS  DURING  USGS- 


Ground  Coverage  During  USGS-3,  Pass  1 
Greensboro  and  Raleigh  1 :250000  Quad  Sheet) 
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Figure  23.  CV-580  SAR  Ground  Coverage  During  USGS-3,  Passes  2  and  3 
(after  USGS  Raleigh  1:250000  Quad  Sheet) 
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3.2  SAR  CALIBRATION  DATA 

The  data  acquired  during  the  US6S-NCDPA  flights  to  assist  in  the 
calibration  of  the  SAR  imagery  are  discussed  in  this  section.  These 
data  include  images  of  a  calibrated  corner  reflector  array  and  images 
internally-generated  calibration  signals.  Appendix  A  to  this  report 
discusses  calibration  of  SAR  imagery  in  detail. 

The  calibrated  corner  reflectors  were  located  in  two  separate 
arrays.  These  arrays  were  located  at  the  Asheboro  Municipal  Airport 
and  at  a  site  in  the  Duke  Forest  east  of  Durham,  North  Carolina. 
The  locations  of  these  two  sites  are  presented  in  Figure  1.  The 
Asheboro  Array  was  imaged  during  three  passes,  and  the  Duke  Forest 
Array  during  six  passes. 

A  total  of  31  precision  corner  reflectors  and  one  active  radar 
calibrator  (see  Brunfelt  and  Ulaby,  1984)  were  deployed  in  the  two 
arrays.  The  types  and  sizes  of  these  reflectors  are  summarized  in 
Table  7,  and  their  placement  within  the  two  arrays  illustrated  in 
Figures  27  and  28.  Figure  29  is  a  photograph  of  typical  reflectors 
used  to  calibrate  SAR  imagery,  and  Figures  30  and  31  contain  X-,  C- 
and  L-band  SAR  images  of  these  two  arrays. 

At  specified  times  during  the  SAR  flights,  internal  calibration 
signals  were  recorded,  which  are  used  to  calibrate  the  SAR  imagery 
(see  discussion  in  Chapter  2).  Figure  32  presents  examples  of  these 
internal  calibration  signals  collected  at  L-band. 

3.3  GEOMETRIC  REFERENCE  REFLECTORS 

One  of  the  goals  of  the  USGS-NCDPA  data  collection  program  was 
to  provide  a  mul ti frequency,  multipolarization,  digitized  SAR  data 
set  to  evaluate  the  utility  of  multi-channel  SAR  data  for  the  detec¬ 
tion  of  stressed  vegetation  communities.  Implicit  in  this  goal  is 
tne  capability  to  merge  the  mul ti -channel  SAR  data. 
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SUMMARY  OF  CALIBRATED  CORNER  REFLECTORS  USED 
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Duke  forest  calibration  array 


■  Active  radar  calibrator 
□  Square  trihedral 
Zl  Triangular  trihedral 


Figure  28.  Approximate  Locations  of  Reference  Reflectors  Within  the  Duke 
Forest  Calibration  Array  (see  Table  7  for  reflector  parameters) 


ALUMINUM  TRIHEDRAL 


STYROFOAM  TRIHEDRAL 


LUNEBERG  LENSE  ACTIVE  RADAR  CALIBRATOR 


Figure  29.  Typical  Calibrated  Reflectors  Deployed  at 
the  Asheboro  and  Duke  Forest  Calibration 
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To  assist  in  this  latter  objective,  a  set  of  specialized  radar 
reflectors  were  designed  and  deployed  around  test  areas  A  and  B. 
These  reflectors,  called  dihedrals  because  of  their  design,  reflect 
the  incident  microwave  radiation  being  transmitted  by  the  SAR  in  such 
a  way  that  a  return  is  detected  by  both  the  like-  and  cross-polarized 
channels  of  the  SAR  system.  The  design  and  construction  of  these 
dihedral  reflectors  is  discussed  in  more  detail  in  Appendix  B,  which 
contains  an  article  which  has  been  submitted  for  journal  publication. 

A  total  of  thirteen  (13)  dihedral  reflectors  were  deployed  at 
nine  different  sites.  Figure  33  presents  a  photograph  of  a  dihedral 
reflector,  and  Figure  34  presents  the  positions  of  the  nine  sites 
where  the  dihedrals  were  deployed.  Table  8  lists  the  latitude  and 
longitude  of  the  positions  of  the  dihedral  reflectors. 

A  review  of  the  SAR  imagery  indicates  that  the  dihedrals  were 
not  detected  every  time  they  were  within  the  SAR  swath.  For  the  X- 
band  data,  the  reflectors  were  detected  72  percent  of  the  time  on 
the  like-polarized  channel  and  81  percent  of  the  time  on  the  cross- 
polarized  channel.  For  the  C-band  data,  the  reflectors  were  detected 
33  percent  of  the  time  on  the  like-polarized  channel  channel  and  56 
percent  of  the  time  on  the  cross-polarized  channel.  Finally,  for 
the  L-band  data,  the  dihedral  reflectors  were  detected  22  percent  of 
the  time  on  the  like-polarized  channel  and  57  percent  of  the  time  on 
the  cross-polarized  channel. 

3.4  DIGITIZED  SAR  DATA 

The  primary  objective  of  this  data  collection  program  was  to  gen¬ 
erate  a  multifrequency  and  multipolarization  SAR  data  set  in  a  digi¬ 
tal  format.  To  produce  this  data  set,  ERIM's  Hybrid  Image  Processing 
Facility  was  used  (see  Ausherman,  et  al.,  1975,  which  is  included  as 
Appendix  C  to  this  report). 


re  33. 
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Scientists  from  the  USGS  selected  imagery  from  twelve  separate 
scenes  to  be  digitized.  The  images  were  digitized  using  an  aperture 
with  an  equivalent  ground  size  of  3  m  by  3  m.  The  data  were  then 
slant-to-ground  range  corrected  and  resampled  to  a  pixel  size  of  3  m 
by  3  m  (see  Section  3.4.1  below).  A  first-order  radiometric  correc¬ 
tion  was  also  applied  to  the  data  (see  Section  3.4.2).  The  original 
uncorrected  data  are  archived  at  ERIM.  The  data  were  then  tran¬ 
scribed  onto  a  set  of  6  CCTs,  with  each  CCT  containing  the  six  chan¬ 
nels  ( X— ,  C-  and  L-band;  VV  and  VH  polarizations)  of  data  from  a 
specific  scene. 

Table  9  summarizes  the  digitized  SAR  data  set.  Presented  in 
Table  9  are  the  CCT  number,  the  file  number,  IPL  tape  number  where 
the  original  data  is  stored,  the  area  of  the  scene,  the  mission  and 
pass  of  the  scene,  the  frequencies  and  polarizations  of  the  SAR  data 
for  that  scene,  the  number  of  lines  per  scene,  and  the  number  of 
pixels  per  line  for  each  scene.  Appendix  D  presents  the  format  of 
the  CCTs. 

As  mentioned  above,  two  corrections  were  applied  to  the  SAR  data 
collected  during  the  USGS-NCDCP.  First,  a  slant-to-ground  range  cor¬ 
rection  algorithm  was  applied  to  the  data.  This  was  followed  by  an 
empirical  radiometric  correction  program  to  remove  the  effects  of 
the  SAR  antenna  pattern,  range  fall-off  and  the  illumination  pattern 
of  the  laser  used  during  digitization  of  the  SAR  data.  Each  of  these 
corrections  are  discussed  iri  a  following  section.  For  an  in-depth 
discussion  of  calibration  of  this  data,  the  reader  is  referred  to 
Appendix  A. 

3.4.1  GEOMETRIC  CORRECTION  OF  SAR  DATA 

Because  of  their  sidelooking  geometry,  SARs  produce  imagery  in 
the  slant  range  plane.  The  across-track  position  of  an  object  in  a 
SAR  image  is  determined  by  its  slant  range,  the  straight  line  dis¬ 
tance  of  the  object  from  the  SAR  platform.  Ground  range,  the 
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distance  of  the  object  from  the  ground  track  of  the  vehicle,  is  the 
representation  desired  by  most  image  interpreters  and  analysts. 

Figure  35  illustrates  the  relationship  between  slant  range  and 
ground  range.  In  actuality,  the  earth's  surface  is  curved,  but  from 
aircraft  altitudes,  the  curvature  is  small  enough  to  be  ignored. 
However,  from  spacecraft  altitudes,  the  earth's  curvature  must  be 
accounted  for. 

In  this  section,  we  will  discuss  the  relationship  between  slant 
and  ground  range  on  a  flat,  horizontal  earth.  The  relationship  be¬ 
tween  slant  and  ground  range  is  non-linear,  but  simple.  The  trigono¬ 
metric  relationship  for  flat  earth  is 

Rg  =  Rs  sin  e  ,  (1 ) 

where  e  is  the  incidence  angle,  Rs  refers  to  slant  range,  and  R^ 
refers  to  ground  range.  The  slant  range  is  always  smaller  than 
ground  range  and  is  highly  dependent  upon  the  angle  at  which  the  ter¬ 
rain  is  viewed.  However,  the  relationship  between  the  change  in 
slant  range  to  the  change  in  ground  range  is 

aRg  =  aRc/sin  e  .  (2) 

The  relationship  between  R  and  R  determines  the  scale  of 
the  image.  Because  this  relationship  differs  at  near  range  compared 
to  far  range,  the  range  scale  of  the  slant  range  image  is  not  con¬ 
stant.  This  type  of  distortion  is  at  its  maximum  when  looking  at 
the  nadir  of  the  vehicle,  where  the  incidence  angle  is  0°.  Near  the 
nadir,  a  comparati vely  small  change  in  slant  range  will  produce  a 
large  change  in  ground  range.  As  the  incidence  angle  is  increased, 
a  smaller  and  smaller  change  will  occur  in  ground  range  for  every 
change  in  slant  range.  At  the  limit  when  the  incidence  angle  is  90° 
(the  radar  beam  is  grazing  the  ground  and  parallel  to  the  ground), 
there  is  no  distortion,  because  the  change  m  slant  range  is  identi¬ 
cal  to  the  change  in  ground  ranee.  Because  of  this  continuous 
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change  in  range  scale  across  the  image,  only  one  range  exists  where 
the  slant  range  scale  equals  the  along-track  scale.  At  nearer  range, 
the  map  scale  will  be  smaller,  and  at  a  farther  range  the  map  scale 
will  be  larger  in  the  slant  range  presentation  of  the  image. 

The  slant-to-ground  range  distortion  can  be  removed  by  digitally 
resampling  the  data  in  the  range  direction.  To  do  this,  the  slant- 
ground  range  relationships  can  be  rewritten  as 

Rg  =  ^R2s  -  h2  ,  (3) 

for  the  flat  earth  or  low-altitude  case,  where  h  is  the  platform 

altitude. 

Because  the  range  and  azimuth  scales  are  not  the  same  after  the 
slant-to-ground  range  conversion  has  been  applied,  the  range  dimen¬ 
sion  is  usually  resampled  to  match  the  azimuth  line  spacing.  Two 

types  of  resampling  techniques  are  commonly  used,  depending  upon  the 

type  of  scene,  the  required  accuracy,  and  the  acceptable  cost  for 

the  correction.  The  simplest  and  least  expensive  resampling  tech¬ 
nique  is  the  nearest-neighbor  method.  In  this  method,  a  grid  is  set 
up  having  equal  intervals  of  ground  range,  aR^,  beginning  at  the 


nearest  range  sample,  R 


That  is,  the  ground  range  at  grid  point 


R  =  R  +  naR  .  (4) 

gn  go  g  v 

The  ground  range  spacing  aR^  is  usually  chosen  to  equal  the  dis¬ 
tance  between  samples  in  the  along-track  direction,  so  that  the  image 
will  have  unity  aspect  ratio  when  displayed  on  a  device  which  has 
pixel  elements  of  equal  size  in  both  dimensions. 

The  resampling  is  then  carried  out  by  considering  each  value  of 
R^n  in  turn,  and  computing  the  corresponding  ground  range  from 
Eq.  (4).  The  input  range  sample  nearest  to  this  computed  ground 
range  is  then  copied  into  the  output  data  set.  The  maximum  error  in 
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the  location  of  any  pixel  using  this  method  is  one  half  of  the  ground 
range  distance  between  pixels.  Usually,  this  is  an  acceptable  error, 
particularly  if  the  scene  reflectivity  distribution  is  relatively 
smooth. 

Another  method  of  resampling  interpolates  between  samples  to  find 
the  data  value  at  the  required  location.  The  same  steps  as  described 
above  are  carried  out  except  that  instead  of  choosing  the  input  range 
sample  nearest  to  the  computed  slant  range  for  the  output  pixel  under 
consideration,  an  interpolation  formula  is  used  which  involves  the 
data  values  at  several  locations  surrounding  this  pixel.  Several 
different  interpolation  formulas  have  been  used,  most  of  them  varia¬ 
tions  or  approximations  of  a  sine  (i.e.,  [sin  x]/x)  function 
(Shuchman,  et  al . ,  1977  ).  This  procedure  results  in  a  theoretically 
more  exact  reconstruction  of  the  image  but  is  more  difficult  to  im¬ 
plement  and  requires  more  computer  time  for  its  operation. 

The  recorded  altitude  (h)  of  the  CV-580  aircraft  was  determined 
by  an  altimeter  which  calculates  altitude  by  measuring  the  difference 
in  barometric  pressure  between  the  aircraft  and  a  point  on  the 
ground,  usually  the  air  field  where  the  airplane  takes  off  and  lands. 
Naturally  occuring  variations  in  atmospheric  pressure  can  result  in 
variations  in  the  actual  altitude  of  the  aircraft.  Studies  have 
shown  that  these  altitude  variations  are  often  as  high  as  50  m 
(Lyden,  1983;  Kasischke,  et  al.,  1983),  which  can  lead  to  errors  in 
ground  range  when  calculated  using  Eq.  (3). 

During  the  North  Carolina  data  collection  flights,  several  cold 
fronts  were  moving  through  the  area,  resulting  in  an  unstable  atmo¬ 
sphere.  It  is  therefore  very  likely  that  the  barometric  pressure 
was  variable,  resulting  in  some  variation  in  the  recorded  aircraft 
alttitude.  The  slant-to-ground  range  corrections  made  on  the  data 
used  the  recorded  altitude,  and  thus  some  ground-range  distortions 
may  be  the  result  of  altitude  variation. 
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If  a  higher  degree  of  accuracy  in  the  slant-to-ground  range  cor¬ 
rection  is  desired,  the  following  method  can  be  used  to  determine 
the  aircraft  altitude  (h). 

If  there  are  two  points  on  a  SAR  image  (see  Figure  36)  whose  lo¬ 
cation  can  be  identified  on  a  map,  then  we  can  easily  calculate  the 

slant  range  distance  to  each  point  (R^  and  R^)  and  the  ground 

range  distance  (in  the  range  direction)  between  the  two  points  (6). 
The  depression  angle  (e^)  to  the  second  point  can  be  calculated 


and  the  altitude  can  be  calculated 


h  =  Rj  sin  eD  .  (6) 

Other  geometric  distortions  will  be  introduced  into  the  image 
because  Eq.  (3)  assumes  a  flat  earth.  Since  there  are  elevation  dif¬ 
ferences  within  the  scene,  these  will  also  lead  to  ground-range  vari¬ 
ations.  These  elevation-induced  variations  can  be  removed  if  the 
SAR  data  are  combined  with  elevation  data  (see  Jackson,  1983). 

3.4.2  RADIOMETRIC  CORRECTION 

A  first-order  radiometric  correction  was  applied  to  the  SAR  data 
in  order  to  produce  an  even  illumination  pattern  across  the  imaged 
swath.  In  this  section,  we  will  first  discuss  the  sources  of  radio- 
metric  variation  in  a  SAR  scene  and  then  discuss  how  these  variations 
were  normalized. 

The  intensity  (Pj)  recorded  on  the  SAR  scene  is  a  function  of 
the  radar  backscatter  (a)  for  the  area  being  illuminated  and  a  vari¬ 
ety  of  radar  parameters.  Assuming  that  the  SAR  is  operating  within 
its  linear  region,  this  relationship  can  be  expressed  as: 
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P-j-G  ( 0 ) X  H^o 

( 4  it  ) 3  R4 


where  Py  is  the  transmitted  power, 

G(e)  is  the  antenna  gain  as  a  function  of  the  incidence  angle, 

e, 

x  is  the  radar  wavelength, 

is  the  radar  system  transfer  function  or  radar  system 
gain, 

R  is  the  slant-range  distance  to  the  target,  and 
Pn  is  the  system  noise. 

Of  interest  in  this  discussion  are  the  effects  of  range  (R)  and 
the  antenna  gain  [G(e)].  If  all  other  SAR  parameters  are  constant, 
a  specific  target  or  feature  will  have  a  lower  P^  the  further  out 
in  slant  range  the  target  is  located  because 

P 1 2 

p —  ~  ~4  <  1  •  (8) 

KI1  R* 

If  we  assume  that  Pp  has  a  value  of  1,  then  the  values  of 

P 1 2 / P 1 1  ^or  increasin9  ranges  are  illustrated  in  Figure  37. 

The  antennas  for  the  X-,  C-  and  L-band  SAR  each  have  a  unique 

pattern,  with  a  peak  gain  or  intensity  occuring  at  the  center  of  the 
antenna  beam.  The  intensity  or  gain  falls  off  as  a  function  of  the 
angular  distance  away  from  the  center  of  the  beam,  i.e.,  a  target's 
intensity  will  be  lower  at  the  outside  of  the  beam's  center  than  it 
will  be  at  the  center  of  the  beam. 

For  the  most  part,  during  the  North  Carolina  data  collection  the 
antennas'  beams  were  centered  at  an  incidence  angle  of  62°. 

Figure  37  also  illustrates  the  expected  ratio  of  PT9/PT1  as  a 


Laser  Illumination  Profile 


Figure  37.  Illustration  of  the  Effects  of  Range  Fall-off,  X-  and 
l-Band  Antenna  Pattern,  and  Laser  Illumination  Pattern 
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function  of  Qj,  where  Pj-|  is  the  intensity  of  a  target  at  the 
center  of  the  antenna  beam.  We  can  see  that  the  L-band  antenna  pat¬ 
tern  is  fairly  flat,  whereas  the  X-band  antenna  fall-off  is  rather 
dramatic.  No  antenna  pattern  is  yet  available  for  the  C-band 
antenna. 

One  more  factor  which  will  affect  the  radiometric  fidelity  of 
the  digital  SAR  data  is  the  illumination  pattern  of  the  laser  beam 
which  was  used  to  optical ly-process  SAR  data  in  the  HIPF.  This  pat¬ 
tern  is  also  shown  in  Figure  37. 

To  be  able  to  quantitatively  compare  the  intensity  values  within 
a  given  scene  or  between  the  same  scenes  of  different  frequencies  or 
between  the  same  area  on  different  passes  requires  calibration  of 
the  data.  This  procedure  is  discussed  in  Appendix  A.  We  have  ap¬ 
plied  a  first  order  radiometric  correction  to  the  data  which  has  es¬ 
sentially  eliminated  any  range  fall-off  which  may  exist  in  the  data. 
This  does  not  imply  the  data  have  been  calibrated. 


The  range  variation  within  a  scene  was  identnied  by  averaging 
approximately  200  pixels  in  azimuth  and  then  plotting  a  running  aver¬ 
age  (of  approximately  30  pixels)  of  pixel  intensity  versus  range. 
This  plot  was  further  smoothed  to  remove  any  obvious  deviations,  and 
a  multiplicative  correction  value,  C-,  calculated  as 


where  i  is  the  pixel  location  in  range,  S.  is  the  plot  intensity 
and  Smax  is  the  maximum  intensity  of  the  curve.  A  corrected  value 
for  each  pixel,  P^,  is  then  calculated  as 

Pii  =  CiPji  •  HO) 

Figure  38  illustrates  an  L-band  image  with  fall-off  in  the  range 
dimension.  Figure  39  shows  the  empirical  radiometric  correction 
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curve  generated  for  these  data.  Finally,  Figure  40  shows  the  radio- 
metrically  and  geometrically  corrected  data. 

3.4.3  DIG  I  TALLY -PROCESSED  SAR  DATA 

Another  method  to  obtain  digital  SAR  image  is  to  digitally  pro¬ 
cess  the  SAR  signal  histories  recorded  on  the  high  density  digital 
tape  (HDDT).  During  the  USGS-NCDPA  program,  two  channels  of  SAR  data 
were  recorded  on  HDDT  during  each  pass  (see  Table  3).  These  data 
are  presently  archived  at  ERIM.  Figures  41  to  43  present  examples 
of  imagery  generated  from  digitally-recorded  data.  These  images  are 
all  from  the  Duke  Forest  area  surrounding  the  calibration  array  (see 
Figure  30) . 


3.5  ENGINEERING  ASSESSMENT 

Immediately  after  optically  processing  the  SAR  imagery  from  a 
data  collection  program,  ERIM  engineers  and  scientists  perform  an 
engineering  assessment  on  the  data.  These  engineers  and  scientists 
have  extensive  experience  in  SAR  image  interpretation  and  recognize 
image  characteristics  which  are  indicative  of  SAR  system  malfunc¬ 
tions.  During  this  engineering  assessment,  a  rigorous  set  of  guide¬ 
lines  is  followed.  The  degree  to  which  each  image  characteri stic  is 
present  is  noted.  The  following  rating  criteria  is  used: 

+  slight  reduction  in  image  quality, 

*  moderate  reduction  in  image  quality,  or 
severe  reduction  in  image  quality. 

In  addition  to  noting  the  types  and  severity  of  the  characteris¬ 
tics  present  in  the  imagery,  the  overall  quality  of  each  pass  of 
imagery  is  rated  using  the  following  categories: 


ERIM 
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7  Excellent  imagery 

6 

5  Good  imagery 

4 

3  Fair  imagery 

2 

1  Poor  imagery 

0  Data  not  recorded  or  of  insufficient  quality  to  be  useful. 

A  rating  of  excellent  implies  that  the  data  are  of  the  highest 
quality  and  require  no  additional  precision  processing.  A  rating  of 
good  implies  that  there  are  slight  or  moderate  image  perturbations 
present,  which  usually  can  be  removed  through  precision  processing 
of  the  imagery.  A  rating  of  fair  implies  there  is  much  useable  data 
present,  and  that  the  moderate  to  severe  image  perturbations  may  be 
removed  through  precision  processing.  A  rating  of  poor  implies  that 
the  severe  image  perturbations  generally  cannot  be  removed  through 
precision  processing. 

Only  two  image  characteristics  were  noted  during  the  engineering 
assessment  of  the  USGS-NCDPA  SAR  data  set,  and  these  were  noted  only 
on  the  X-band  data.  These  two  characteristics  were  modulation  inten¬ 
sity  variations  and  film  transpart  runaways. 

Excessive  aircraft  motion  caused  by  areas  of  high  turbulence 
during  data  collection  sometimes  create  a  range-oriented  banded  pat¬ 
tern  within  the  imagery.  These  bands  are  referred  to  as  modulation 
intensity  variations.  This  type  of  image  character i Stic  is  most 
often  found  in  X-band  data.  Figure  44  presents  an  example  of  a  set 
of  modulation  intensity  variations  from  USGS-2,  Pass  7. 

A  periodic  malfunction  of  the  camera  which  recorded  X-band  data 
resulted  in  gaps  in  the  recorded  SAR  signal  histories.  These  gaps 
result  in  blank  areas  on  the  output  imagery.  Figure  44  presents  an 
example  of  X-band  imagery  where  this  film-drive  runaway  occurred 


(a)  Modulation  Intensity  Variation 
(USGS-2,  Pass  7) 


(b)  Film  Transport  Runaway 
(USGS-2,  Pass  7) 


Figure  44.  Examples  of  Image  Characteristics  Found  on  USGS-NCDPA 
X-Band  SAR  Data 
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during  USGS-1,  Pass  7.  No  data  are  lost  due  to  these  runaways,  but 
the  optically  processed  imagery  from  passes  where  this  problem  oc¬ 
curred  will  have  to  be  mosaiced  in  order  to  produce  a  continuous 
image. 

The  image  characteristics  noted  on  the  X-band  SAR  imagery  were 
all  slight  in  nature,  and  did  not  significantly  reduce  the  utility 
of  the  data  set.  Overall,  the  SAR  data  set  collected  over  the  NCDPA 
was  of  extremely  high  quality.  The  average  quality  rating  for  the 
X-band  data  was  6.2,  for  the  L-band  data,  7.0,  and  for  the  C-band 
data,  7.0.  No  additional  optical  processing  is  required  for  this 
data  set. 


SUMMARY 


8etw een  8  and  10  April,  1984,  three  SAR  missions  were  flown  over 
the  North  Carolina  Digital  Project  Area  (NCDPA)  by  the  ERIM/CCRS 
CV-580  SAR  System  in  order  to  collect  mu  1 ti frequency  and  multi- 
polarization  SAR  imagery  over  selected  USGS  test  sites.  A  total  of 
22  SAR  passes  were  collected  over  five  different  test  sites.  Approx¬ 
imately  790  line  kilometers  (7900  square  kilometers)  of  high  quality 
SAR  data  were  collected. 

These  data  were  first  optically  processed  into  imagery.  Copies 
were  provided  to  USGS  personnel,  who  then  selected  twelve  areas  to 
be  digitized.  A  digitized  SAR  data  set  consisting  of  three  frequen¬ 
cies  ( X— ,  C-  and  L-bands)  and  two  polarizations  for  each  frequency 
(VV  and  VH)  was  then  produced  on  ERIM's  Hybrid  Image  Processing  Fa¬ 
cility.  First-order  geometric  and  radiometric  corrections  were  ap¬ 
plied  to  the  data.  Thus,  a  data  set  has  been  generated  to  use  in 
the  evaluation  of  SAR  imagery  for  detection  of  geologically-stressed 
vegetated  areas. 

In  order  to  realize  the  full  potential  of  the  NCDPA  SAR  data  set, 
users  should  realize  that  sophisticated  analysis  may  need  to  be  em¬ 
ployed.  First,  to  draw  quantitative  conclusions  from  the  changes  in 
radar  backscatter  from  different  passes,  the  SAR  data  must  first  be 
calibrated.  Calibration  of  the  SAR  data  requires  reduction  of  the 
calibration  array  data  and  calibration  signal  data  as  well  as  fol¬ 
lowing  a  careful  set  of  procedures  (see  Appendix  A  or  Larson,  et  al., 
1985). 

Second,  higher-order  geometric  corrections  (see  Jackson,  1983) 
may  be  applied  to  the  original  SAR  data  in  order  to  reduce  radar  lay¬ 
over  and  radar  foreshortening.  These  higher-order  geometric  correc¬ 
tions  would  not  only  make  the  data  easier  to  merge,  but  also  could 
be  used  to  account  for  scattering  variations  induced  by  changes  in 
the  local  angle  of  incidence. 
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Finally,  if  analysis  of  the  hybrid-digital  data  shows  significant 
results,  researchers  should  consider  the  use  of  digitally-processed 
SAR  data.  These  data  have  a  higher  spatial  resolution,  lower  noise 
levels,  and  higher  dynamic  range  than  the  hybrid-digital  data.  In 
addition,  image  speckle  can  be  reduced  through  multiple-look 


processing. 
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APPENDIX  A 

CALIBRATION  OF  SAR  DATA 


Calibration  of  SAR  data  is  a  relatively  strai ght-forward  concept 
to  understand  in  a  theoretical  sense,  although  the  implementation  of 
a  total  calibration  algorithm  requires  a  careful,  routinized  proce¬ 
dure  in  order  to  keep  track  of  all  the  factors  which  can  affect  each 
step  in  the  algorithm.  In  this  Appendix,  we  will  discuss  the  funda¬ 
mentals  of  SAR  image  calibration.  For  a  more  in-depth  discussion  of 
calibration,  the  reader  is  referred  to  Walker  and  Larson  (1981), 
Larson,  et  al.  (1 982,  1985)  and  Larson  and  Maffett  (1985).  In  our 
discussion,  we  will  use  data  collected  during  a  SAR/ocean  imaging 
experiment  conducted  for  the  Office  of  Naval  Research  (see  Kasischke, 
et  al.,  1985)  to  illustrate  the  basic  concepts  behind  SAR  image  cali¬ 
bration.  Figure  A-l  presents  X-  and  L-band  SAR  data  of  internal  wave 
surface  patterns. 

The  goal  of  the  calibration  algorithm  is  to  relate  the  measured 
intensity,  Pj,  from  a  processed  SAR  image  to  radar  cross-section 
(a).  Assuming  the  SAR  is  operating  within  its  linear  region,  this 
relationship  can  be  expressed  as 

PTG2(e)x2H  a 

P^- - 3“4^~  +  pn  (A-D 

1  (4*)  R 


where 


Pj  is  the  transmitted  power, 

G(e)  is  the  antenna  gain  as  a  function  of  the  incidence  angle, 


x  is  the  radar  wavelength 

is  the  radar  system  transfer  function  or  radar  system 


R  is  the  slant-range  distance  to  the  target,  and 
3  is  the  system  noise. 
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Using  SAR  imagery  of  calibrated  reflectors  of  known  radar  cross- 
section  (o)  and  internal  calibration  signals,  the  basic  calibration 
approach  used  at  ERIM  is  as  follows: 

1.  Inspect  the  candidate  SAR  imagery  for  artifacts  and  satura¬ 
tion  to  determine  its  potential  for  calibration, 

2.  Digitally  process  the  SAR  data  set, 

3.  Calculate  H  and  P„  for  the  scene  to  be  calibrated  using 

s  n  3 

the  calibration  signals  for  that  scene  and  the  calibration 
reflector  scene,  and 

4.  Calculate  a  for  the  image  scene  using  P^  from  that  scene 

and  the  specific  parameters  for  that  scene  (e.g.,  a  new  Pj 
and  G(e)  and  R). 

Figure  A-2  outlines  the  overall  procedure  used  to  calibrate  SAR 
imagery.  The  first  two  steps  (inspection  and  processing  of  the  SAR 
data)  appear  to  be  straightforward,  but  in  fact  are  crucial  steps  in 
the  analysis  and  calibration  process.  The  importance  of  this  step 
cannot  be  overemphasized. 

In  selection  of  a  SAR  data  set  to  address  a  specific  analysis 

topic,  it  is  necessary  to  determine  if  the  SAR  signatures  are  within 
the  linear  response  region  of  the  SAR  system  as  well  as  determine 

whether  or  not  the  calibration  signals  are  adequate  for  use  in  deter¬ 
mining  the  SAR  system  gain.  Once  these  two  questions  have  been  ad¬ 
dressed,  the  SAR  data  can  be  processed  into  imagery  and  the  calibra¬ 
tion  procedure  begun. 

Figure  A- 3  presents  image  intensity  scans  through  the  internal 
wave  surface  patterns  present  in  Figure  A-2.  The  scan  was  made  in 
the  range  direction,  with  the  ship  present  in  the  image  being  the 
center  of  the  scan.  The  scan  was  generated  by  averaging  100  lines 
in  azimuth  and  20  elements  in  range.  It  is  these  intensity  values 
which  will  be  used  to  illustrate  the  calibration  process. 
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igure  A-3.  Uncal i brated  SAR  Image  Intensity  Scans  Through 

Surface  Internal  Wave  Patterns  From  SARSEX-4,  Pass  5 


An  examination  of  the  intensities  of  the  SAR  signal  histories 
will  indicate  whether  the  SAR  system  is  operating  in  the  linear  re¬ 
gion  in  the  area  of  interest,  i.e.,  whether  the  signals  are  saturated 
or  not.  To  do  so  requires  comparing  the  signals  from  an  area  of 
interest  to  a  curve  of  the  intensities  from  the  calibration  signals, 
some  of  which  are  saturated.  The  intensities  of  the  SAR  signal  his¬ 
tories  are  illustrated  in  Figure  A-4,  which  presents  a  curve  gener¬ 
ated  from  the  internal  calibration  signals  along  with  the  location 
of  the  intensities  of  the  ship,  background  level  and  brightest  inter¬ 
nal  wave  value. 

The  second  step  in  the  calibration  procedure  is  to  calculate  each 
of  the  factors  in  Eq.  (A-l).  This  includes  using  the  SAR  imagery  of 
the  set  of  calibrated  reflectors.  Once  the  SAR  images  of  the  re¬ 
flectors  are  processed,  an  intensity,  P,  is  extracted  for  each  corner 
reflector  in  the  scene.  This  intensity  is  calculated  as 

n  m 

p  •  E  Z  («> 

1=1  j=l 

where  x.^.  is  the  power  of  a  pixel, 
i  is  the  line  number,  and 
j  is  the  element  number. 

It  is  then  necessary  to  calculate  the  background  clutter  for  the 
area  where  the  reflectors  have  been  deployed.  Once  this  clutter 
value  (C)  has  been  calculated  and  normalized  to  the  same  area  as  the 
area  used  for  the  reflector,  then  Pj  is  calculated 

Pj  =  P  -  C  ( A-3 ) 

It  should  be  noted  that  the  same  number  of  elements  (i)  and  lines 
(j)  should  be  used  to  calculate  P^  for  each  reflector.  Figure  A-5 
presents  the  P^  vs.  a  curves  generated  for  the  calibration  array 
imaged  during  the  ONR  experiment. 
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Through  the  establishment  of  a  linear  relationship  between  o  and 
Pj,  we  can  now  arrange  Eq.  (A-l)  to  calculate  Hg  for  the  cali¬ 
brated  reflector  array.  Specifically 


P  T ( 4w ) 3R4 

H  =  — ~ r 
s  PTG^(e)x  o 


( A— 4 ) 


Table  A-l  lists  the  various  factors  to  calculate  H$  for  the 
calibration  array. 

The  next  step  in  the  calibration  algorithm  is  to  calculate  the 
system  gain,  H  ,  for  the  imaged  scene.  This  is  accomplished 
through  measurement  of  the  calibration  signals  generated  before  or 
after  each  pass  of  SAR  imagery.  These  signals  are  measured  in  much 
the  same  manner  as  the  reflector  images  using  Eq.  (A-3).  However, 
instead  of  measuring  a  clutter  factor,  a  system  noise  factor  is  cal¬ 
culated  and  subtracted  from  the  power  value  (see  Lyzenga  and 
Shuchman,  1984).  These  power  values  are  then  plotted  against  the 
input  power  used  to  generate  the  signal.  A  curve  is  generated  using 
both  sets  of  calibration  signals  (for  the  calibration  array  and  the 
SAR  image  to  be  calibrated).  The  difference  between  the  linear  por¬ 
tions  of  these  curves  (k)  is  the  system  gain  difference  between  the 
two  passes.  The  calibration  signal  curve  at  X-band  for  the  two  ONR 
data  sets  are  presented  in  Figure  A-6.  The  system  gain  for  the  image 
scene,  H$  calculated  from  Eq.  (A-4),  is  thus  calculated  by  adding 
or  subtracting  k  from  H$.  For  our  example,  k  equals  18.0  dB  at 
L-band  and  18.5  dB  at  X-band. 

Finally,  a  radar  cross  section  for  any  pixel  in  the  image  scene 
can  be  calculated 


(Pj  -  Pn)(4n)3R4 

PTG2(e)xV 


( A— 5 ) 
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TABLE  A-l 

SYSTEM  PARAMETERS  USED 
SAR  IMAGERY 

TO  CALIBRATE 

a. )  Parameters  for 
from  SARSEX  4, 

Calculation  of  Hs  for 
Pass  12 

Calibration  Reflectors 

Parameter 

X-Band 

L-Band 

PI  -Pn 

0 

76,504 

5,248 

X 

0.032  m 

0.24  m 

R 

9,666  m 

9,516  m 

G2(e) 
for  e  = 

48.6  dB 

28.2  dB 

PT 

1,741  W 

5,650  W 

HS 

1.0301  x  1019  W-1 

3.9974  x  1017  W"1 

b.)  Parameters 

for 

calculation  of  o  values 

for  SARSEX-4,  Pass  5 

P 

n 

11,442 

204 

X 

0.032  m 

0.24  m 

R 

7,675  to  1 1 ,785  m 

7,675  to  11,785  m 

9I 

29°  to  55° 

29°  to  55° 

9j (Peak ) 

62° 

62° 

G2(e) 

Variable  Dependent  on 

9j  and  9 j (Peak ) 

PT 

1,700  W 

5,452  W 

H 

7.2923  x  1020  W-1 

19  -1 

2.5222  x  10  y  M 

•ltVW 
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where  Py  if  the  transmitted  power  and  Pn  is  the  system  noise  for 
that  scene,  and  R  and  G(e)  are  calculated  specifically  for  that 
pixel . 

Using  Eq.  (A-5),  a  radar  cross-section  can  be  calculated  for  any 
pixel  or  area  of  the  imaged  scene.  Figure  A-7  presents  the  cali¬ 
brated  radar  cross-section  curves  for  the  ONR  SAR  imagery. 
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Abstract  -  Dihedral  reflectors  can  be  used  for  both  cross-  and 
like-polarized  radar  backscatter  to  calibrate  intensity  and  provide 
locational  accuracy  for  image  regi strati  on .  They  can  be  economically 
constructed  so  that  a  large  number  can  be  used  for  calibration  pur¬ 
poses  across  an  imaging  radar  swath.  Both  rotating  platform  and  SAR 
data  demonstrated  the  dihedral's  partitioning  of  reflected  energy 
between  cross-  and  like-polarization. 


1.  INTRODUCTION 

Research  in  the  uses  of  synthetic  aperture  radar  (SAR)  data  is 
focusing  more  and  more  on  measuring  absolute  and  relative  radar  back¬ 
scatter  as  well  as  merging  mu  1 ti -frequency  and  mul ti -pol ar ization 
SAR  data  collected  during  different  SAR  overflights.  Calibration  of 
SAR  data  requires  the  use  of  calibrated  reference  reflectors  and  in¬ 
ternal  calibration  signals  [1].  Merging  of  multiple  SAR  data  sets 
is  aided  by  using  reference  reflectors  whose  positions  are  known. 

Most  SAR  calibration  experiments  to  date  have  used  precision 
trihedral  (corner)  reflectors  whose  cross-section  has  either  been 
calculated  or  measured.  Recently,  Brunfeldt  and  'Jlaby  nave  developed 
an  active  raoar  calibrator  (ARC)  which  is  essentially  a  transponder 
[ZJ.  The  ARC  can  be  maoe  to  generate  either  like-  or  cross-polarized 
returns  or  both  simultaneously. 
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An  increasing  amount  of  research  is  now  being  conducted  on  the 
analysis  of  cross-polarized  SAR  imagery.  Trihedral  reflectors  do 
not  return  cross-polarized  waves  and  thus  cannot  be  used  as  a  cali¬ 
bration  source.  ARCs  can  be  used,  but  since  it  is  necessary  to  have 
a  large  number  of  different  sized  calibration  reflectors  to  generate 
a  proper  calibration  curve,  it  would  be  quite  costly  to  deploy  a  suf¬ 
ficient  number  for  calibration  purposes. 

We  wish  to  describe  a  dihedral  reflector  which  returns  both 
like-  and  cross-polarized  EM  waves.  The  dihedral  reflector  is  rela¬ 
tively  easy  to  construct  and  deploy,  and  cost-effective  for  a  large 
number  of  emplacements.  This  type  of  reflector  can  be  used  to  cali¬ 
brate  both  like-  and  cross-polarized  SAR  imagery  as  well  as  to  serve 
as  a  position  reference  to  merge  multiple-date  SAR  data  sets. 

II.  DESIGN  AND  EMPLACEMENT 

A  polarized  wave  striking  a  conducting  surface  at  a  given  angle 
<t>  will,  upon  reflection,  be  rotated  in  polarity  by  this  angle.  Two 
successive  reflections,  such  as  occur  from  a  dihedral  reflector,  will 
cause  the  polarity  to  rotate  by  2 i>.  If  the  dihedral  is  oriented  so 
that  the  intersection  of  the  plates  is  0  =  22  1/2°  with  respect  to 
the  direction  of  the  polarity  and  is  perpendicular  to  the  propagation 
direction  of  an  incident  radar  beam,  the  polarity  of  the  reflected 
beam  will  be  rotated  by  45°.  Since  a  45°  rotation  has  equal  compo¬ 
nents  of  like-  and  cross-polarizations,  both  polarizations  can  be 
sensed  in  the  radar  return. 

Figure  1  schematically  illustrates  a  dihedral  reflector,  which 
consists  of  two  joined,  perpendicular  sheets  of  reflecting  material 
with  dimensions  a  and  b.  The  proper  orientation  of  the  dihedral  with 
respect  to  the  polarity  of  the  radar  wave  can  be  maintained  by  ad¬ 
justing  the  heights  of  its  front  and  back  corners  above  a  flat  sur¬ 
face.  These  heights  (see  Figure  1)  are  adjusted  as  follows: 


99 


A  =  (a//2)  sin  e  , 


and 


8  =  (ih'7)  cos  9  sin  it  , 

where  a  is  the  dimension  of  the  dihedral  plate  as  shown  in  Figure  1, 
9  is  the  radar  depression  angle,  and  it  is  half  the  angle  of  the  de¬ 
sired  polarity  rotation  upon  reflection. 


Then,  for  e  >  ^ 

h1  =  A  +  8  , 
h2  =  A  -  B  , 


and  for  9  <  it 


h1  =  2A  , 
h2  =  0  , 


where  the  heights  h^,  h^,  and  hfa  are  shown  in  Figure  1.  The 
dimensions  h^  and  h^  can  be  interchanged  with  the  amount  of  rota¬ 
tion  remaining  constant  but  the  direction  being  reversed.  The  back- 
scatter  cross-section  of  such  a  dihedral  [3]  is 


_  8naV  , 
°max  “  ^2 


where  r  is  the  maximum  cross-sect  ion,  and  \  is  the  radar 
ma  x 

wavelength . 


The  dihedrals  can  be  constructed  of  sheet  aluminum  attached  to  a 

wooden  frame,  as  shown  in  Figure  2.  This  dihedral  (3  ft  (0.9  m)  x  4 

2  2 

ft  (1.2  m) )  has  a  of  566  m  for  L-band  and  34,700  m  for 

X-band.  It  is  readily  emplaced  with  proper  orientation  in  the  field 
by  tying  the  corners  to  construction  stakes.  To  aid  in  field  em¬ 
placement  the  proper  heights  of  the  dihedral  corners  can  be  listed 
in  a  table  of  SAR  depression  angles. 

Prior  to  deployment  of  the  dihedral  reflector  shown  in  Figure  2, 
laboratory  measurements  were  obtained.  A  rotating  platform  in  the 
field  of  a  fixed  antenna  was  used  for  this  test  [4].  Figure  3  illus¬ 
trates  pol arity-parti tioned  radar  returns  from  a  dihedral  reflector 
oriented  at  22.5°  with  respect  to  the  incident  polarity  when  the 
dihedral  faces  the  transmitter.  The  dihedral  was  rotated  in  the 
field  of  a  transmit-recei ve  antenna,  and  like-  and  cross-polarized 
radar  returns  plotted.  The  angle-dependent  returns  from  both  polari¬ 
zations  are  equal  or  very  similar  in  value  within  a  rotating  platform 
angle  of  7°  on  either  side  of  the  maximum  return. 

III.  FIELD  TEST 

A  set  of  dihedral  reflectors  was  deployed  during  a  SAR  data  col¬ 
lection  program  conducted  for  the  U.S.  Geological  survey  in  north- 
central  North  Carolina  in  April  1984.  The  ERIM-CCRS  X-C-L  SAR  [5] 
was  used  for  this  purpose.  X-band  (3.2  cm)  and  L-band  (23.5  cm)  data 
collected  over  the  Duke  University  Forest  near  Durham,  North 
Carolina,  demonstrated  the  effectiveness  of  these  dihedrals  for  sens¬ 
ing  both  like-  and  cross-polarized  radar  backscatter.  Figure  4  shows 
a  portion  of  this  imagery  in  both  like-  and  cross-polarized  SAR 
images.  The  dihedral  reflector  shown  in  Figure  2  was  located  near 
the  center  of  each  image,  and  is  evident  in  both  like-  and  cross- 
polarized  images  at  both  wavelengths. 

Note  that  the  'mage  of  the  dihedral  is  stronger  on  the  cross- 
polarized  than  on  the  1 i ke-pol azar i zed  image,  indicating  that 


the  dihedral  reflector  was  deployed  at  an  angle  greater  than  22.5° 

with  respect  to  the  incident  polarity  direction  at  the  actual  SAR 

2 

depression  angle.  Note  also  that  the  large  (8300  m  )  trihedral 
reflectors,  which  return  only  like-polarized  waves,  are  just  visible 
on  the  X-band  cross-polarized  image.  These  reflectors  produced  such 
large  like-polarized  returns  that  they  were  detected  across  the  -23 
dB  isolation  between  the  like-  and  cross-polarized  channels. 
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FIGURES 


Schematic  illustration  of  a  dihedral  reflector. 

Field  emplacement  of  dihedral  reflector. 

Like-Polarized  (Solid  Line)  and  cross-polarized  (Dashed  Line) 
backscatter  vs  angle  of  incidence.  Rotating  platform  in  field 
of  fixed  X-band  antenna.  Aluminum  dihedral  (a  =  6  in,  b  *  12 
in)  positioned  for  22.5°  angle  between  incident  angle  and  line 
of  dihedral  plate  intersection. 

L-Band  and  X-Band  SAR  data  collected  over  Duke  University  Forest 
in  April,  1984.  Dihedral  reflector  is  evident  on  all  four  like- 
and  cross-polarized  images.  Small  corner  reflectors,  located 
above  and  below  the  dihedral  reflectors  are  evident  on  the  like- 
polarized  X-band  image. 
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ABSTRACT 

A  system  is  described  which  is  dedicated  to  Che 
signal  processing  and  image  exploitation  aspects  of 
synchecic  aperture  radar  (SAR).  The  sys tem ut ilizes 
coherent  optical,  hybrid  opt ical-d igital ,  and  wholly 
digital  approaches  to  SAR  data  processing.  The  sub- 
svstems  using  each  one  of  these  processing  media 
are  brieflv  reported.  Examples  of  processed  ima¬ 
gery  are  included  and  an  application  to  radar  remote 
sensing  of  the  earth's  surface  is  given  as  an 
example  of  one  use  of  the  facility. 


INTRODUCTION 

A  svstem  has  been  developed  which  is  dedicated 
to  the  signal  processing  and  image  exploitation 
aspects  of  svnthetic  aDerCure  radar  (SARI.  The 
overall  system  is  capable  of  utilizing  coherent  op¬ 
tical,  hvbrid  opt  ion! -digital .  and  wholly  digital 
processing  technologies  in  the  generation,  enhance¬ 
ment,  and  subsequent  display  of  radar  imagery,  ihe 
ERIM  Radar  Data  Processing  and  Exploitation  Facil- 
itv  has  proven  applicability  in  many  imaging  radar 
app 1 ica t ions ,  including  the  remote  sensing  of 
earth's  resources. 

The  concept  of  Svnthetic  Aperture  Radar  is  well 
documented! 1-3] .  The  most  common  SAR  utilizes  the 
sidelooking  radar  configuration  shown  in  Figure  1. 
The  objective  of  such  a  system  is  to  generate  a 
two-dimensional  continuous  image  of  the  terrain 
strip  illuminated  bv  a  microwave  beam  as  the  radar 
vehicle  moves  in  the  azimuth  I along- t rack)  direc¬ 
tion  at  a  constant  velocity.  Range  ( cross- t rac k ) 
resolution  is  obtained  by  transmitting  disperses 
pulses  and  applying  pulse  compression  techniques 
to  the  returned  signals.  Azimuth  resolution  is 
obtained  bv  recording  the  Doppler  rrecuencv  shi: ts 
associated  with  the  returns  from  poinr  stutterers 
as  then  migrate  through  the  antenna  beam.  Know¬ 
ledge  of  the  Doppler  freauencv  vs.  time  relation¬ 
ship  for  a  point  scatterer  at  a  given  ranee  allows 
one  to  precisely  locate  the  scatterer  m  azimuth 
in  a  manner  analogous  to  the  tmlse  compression 
applied  in  the  range  direction.  The  uncompressed 
radar  echoes  are  usually  recorded  ana  subsequent 
two-dimensional  signal  processing  is  applied  to 
form  an  image.  The  signal  storage  medium  is  usua.iy 
optical  analog  (on  film)  or  digital,  depending  on 
the  processing  medium  Co  be  used  tor  sicr.u.  im¬ 
pression.  The  svstem  described  herein  ilso  i._  ws 


optical  recording  and  subsequent  digital  processing 
of  the  signal  histories. 

For  coherent  optical  signal  compression,  film 
is  che  usual  storage  medium.  The  radar  returns  are 
coherently  mixed  to  video  frequencies  and  used  Co 
modulate  the  intensity  of  a  CRT  spot  as  its  image 
is  swept  across  a  photographic  film  at  a  conscanc 
velocity.  The  video  signal  maintains  a  small  car¬ 
rier  frequency  to  avoid  any  spectral  folding,  and 
it  is  recorded  about  a  D.C.  bias  level  to  enable 
the  proper  recording  of  bipolar  signals.  The  re¬ 
sultant  film  recording  is  a  linear  correspondence 
between  slant-range  (time-delay)  and  the  spatial 
dimension  across  the  film.  Subsequent  pulses  are 
recorded  side-by-side  along  che  film  as  the  film 
is  moved  past  che  CRT  at  a  rate  proportional  to  the 
radar  vehicle  velocity.  The  spatial  dimension 
along  the  signal  film  is  therefore  linearly  related 
to  the  azimuth  dimension  of  the  terrain.  The  de¬ 
veloped  film  forms  the  input  co  che  optical  pro¬ 
cessor. 

OPTICAL  RADAR  SICNAL  PROCESSING 

The  optical  processor  acts  as  a  two-dimensional 
filter  matched  to  the  signal  generated  by  a  point 


FIGURE  1.  Schematic  Representation  o:  Strip- .'.a? 
Mode  Imagine  Radar. 
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FIGURE  2.  Tilted  Plane  Optical  Processor 


scatterer.  Thus,  the  required  tuo-d imens Ional  pulse 
compression  can  be  viewed  as  the  convolution  of  an 
appropriate  reference  function  with  the  recorded 
radar  signal.  In  the  optical  processor  this  convo¬ 
lution  is  performed  by  configuring  che  lenses  such 
that  the  impulse  response  of  the  processor  is  iden¬ 
tical  to  the  desired  reference  function.  The  pri¬ 
mary  advantage  of  the  optical  approach  is  the  nearly 
instantaneous  processing  capability  for  two-dimen¬ 
sional  signals  with  very  high  information  content. 

The  optical  processing  facility  usually  utilizes 
a  configuration  known  as  the  tilted  plane  optical 
processor.  This  optical  system  has  been  previously- 
reported  [1]  and  will  not  be  described  in  detail 
herein.  For  a  brief  description  consider  the  sche¬ 
matic  diagram  of  the  processor  shown  in  Figure  2. 

A  spherical  telescope,  usually  of  unitv  magnifica¬ 
tion,  is  comprised  of  spherical  lenses  and  L-,. 

A  cylindrical  celescope  of  variable  magnification, 
consisting  of  lenses  and  L,  Is  in  tandem.  The 
signal-film  transparency  S  is  illuminated  with  col¬ 
limated,  nonochromatic  (laser)  light  by  utilization 
of  a  point  source  P  and  collimator  Ln-  The  desired 
image  forms  at  the  output  plane  I,  where  it  can  be 
viewed  directly  or  recorded  on  photographic  film. 
Frequency  plane  filtering  can  easily  be  accomplished 
by  placing  che  appropriate  slits  and  weighting  masks 
in  plane  F,  which  is  the  spatial  two-dimensional 
Fourier  transform  of  che  input  data.  Such  filtering 
is  used  primarily  to  limit  the  passbana  of  the  pro¬ 
cessor.  Fcr  certain  forms  of  coded  radar  pulses 
the  desired  matched  filtering  occurs  as  a  conse¬ 
quence  of  the  lens  action  and  free  space  transforma¬ 
tions  associated  with  the  optical  configuration. 

If  che  ranging  signal  transmitted  by  the  radar 
is  a  linear  frequency  modulated  (chirp)  pulse,  the 
slant-range  image  is  formed  near  the  signal  film 
due  to  the  propagation  effects  between  the  film  and 
the  image.  This  is  analogous  Co  the  formation  of  a 
h-  lographic  image.  This  image  is  then  re-imaged  via 
the  spneric.il  celescope  to  the  output  plane  I.  The 
azimuth  signal  history  is  a  true  microwave  holo¬ 
gram  [ 3  which  forms  an  image,  usually  in  a  differ¬ 
ent  plane  and  of  different  magnif ication  than  the 
range  image.  3y  adjusting  che  magn i f  ig at  ion  and  the 
position  of  the  cylindrical  Celescope  (with  no 
effect  on  the  rir.ge  image:,  the  range  and  azimuth 
image  planes  :  an,  be  brought  into  coincidence  with, 
unity  aspect  at  plane  I.  For  processing  of  long 
image  strips  both  input  and  output  films  must  be 


translated  through  che  processor  in  synchronism, 
forming  a  continuous  convolution  in  the  azimuth  di¬ 
rection. 

The  above  processing  method  is  complicated  in 
that  the  azimuth  image  plane  is  tilted  with  respect 
to  the  signal  film  while  the  range  image  is  parallel 
to  the  signal  film1'.  The  solution  is  to  tilt  the 
input  and  output  planes  the  proper  amount  to  bring 
the  planes  to  a  common  focus  [4] .  When  the  input 
plane  is  rotated  abouc  the  azimuth  coordinate,  the 
slant  range  image  formed  by  the  unity  telescope  ro¬ 
tates  the  same  amount.  The  azimuth  image,  however, 
is  demagnified  by  che  cylindrical  telescope  and 
therefore  rotates  a  proportionally  lesser  amount.  A 
tilt  angle  usually  exiscs  which  brings  the  two  image 
planes  into  focus. 

Frequency  plane  filtering  is  required  in  the 
processor  to  remove  the  conjugate  (the  negative  fre¬ 
quency  band)  and  undiffracted  (D.C.)  wavefronts 
which  result  from  the  holographic  signal  recording. 

A  frequency  plane  slit  is  used  to  insure  that  these 
signal  components  do  not  corrupt  the  final  image. 

In  addition,  aperture  weighting  masks  can  be  in¬ 
serted  in  the  frequency  plane  to  provide  the  desired 
processor  impulse  response. 

The  use  of  film  as  the  hardcopy  image  storage 
medium  limits  the  oucpuc  dynamic  range  of  che  system 
to  approximately  25  dB.  The  next  section  will  de¬ 
scribe  a  hybrid  optical-digital  system  capable  of 
providing  direct  digitization  of  the  output  image 
with  increased  dynamic  range. 

HYBRID  OPTICAL-DIGITAL  RADAR  SIGNAL  PROCESSING 

The  image  dissector  and  digitization  (IDD>  faci¬ 
lity  consists  of  a  coherent  optical  processor  inter¬ 
faced  to  a  computer  system  which  provides  digitiza¬ 
tion,  recording,  and  iisplav  of  processor  outputs. 
The  hybrid  concept  effectively  combines  the  power 
and  speed  of  coherent  optical  data  processing  with 
the  flexibility  and  precision  of  digital  processing 
methods.  The  system  can  be  used  for  digitization 


* 

The  azimuth  image  plane  is  tilted  due  to  the 
fact  that  the  Doppler  frequency  vs.  azimuth  posi¬ 
tion  relationship  is  a  function  of  the  range 
the  scatterer.  This  in  turn  implies  that  the  azi¬ 
muth  holographic  signals  have  focal  lengtr.s  which 
vary  linearly  with  range. 
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of  film-recorded  SAR  signal  histories  or  of  opti¬ 
cally  processed  imagery  having  a  dynamic  range  of 
up  to  40  dB.  If  the  digital  phase  histories  are 
digitally  processed  inco  imagery,  it  can  result  in 
a  very  large  dynamic  range  (50  to  60  dB)  previously 
only  available  as  a  light  distribution  in  the  output 
plane  of  an  optical  processor.  Such  imagery  is 
ideal  for  subsequent  exploitation  processing.  The 
unique  design  parameters  for  this  system  were  dyna¬ 
mic  range  response,  resolution,  linearity,  and  total 
system  stability. 

A  block  diagram  of  the  hybrid  facility  is  shown 
in  Figure  3.  The  key  to  the  optical-digital  inter¬ 
face  is  Che  image  dissector  camera  manufactured  by 
ITT  Aerospace/Optical  Division.  The  image  dissec¬ 
tor  has  been  described  as  a  photomultiplier  with  a 
small  electronically  movable  phococathode  area,  thus 
acting  as  an  all-electronic,  low-inertia  microphoco- 
meter.  The  dissector  camera  is  positioned  such  that 
the  oucput  image  of  the  opcical  processor  is  focused 
upon  the  phococathode.  Electrons  are  emitted  from 
the  back  of  the  photocathode,  forming  an  electronic 
image  with  current  density  modulated  according  to 
the  image  input.  The  electron  image,  focused  by 
magnetic  deflection,  falls  upon  an  aperture  plane 
at  the  ocher  end  of  the  drift  tube.  The  aperture 
samples  the  image  by  allowing  only  a  small,  well 
defined  area  of  the  electron  image  to  pass  through 
to  the  electron  multiplier,  which  Chen  multiplies 
the  sampled  pjaocoelectrons  by  a  factor  of  approxi¬ 
mately  5  x  10  .  The  entire  electron  image  is  de¬ 
flected,  allowing  the  aperture  to  sample  different 
points  in  the  input  image. 

In  Che  IDD  system  the  aperture  is  "scanned"  only 
in  the  range  direction.  The  aperture  is  effectively 
dispiaceu  in  che  azimuth  direction  as  the  indue  film 
is  translated  between  successive  range  scans.  The 
dissector  aperture  used  is  rectangular,  measuring 
0.01  mm  in  the  range  direction  and  0.5  mn  in  the 
azimuth  direction.  The  effective  width  of  che 
azimuth  aperture  can  be  reduced  by  placing  a  slit  of 
the  desired  width  in  the  output  plane  of  the  optical 
processor  and  then  re-imaging  the  output  on  the  dis¬ 
sector  phococathode  surface.  The  spatial  frequency 
response  of  che  dissector  exhibits  approximately  40% 
modulation  at  30  line  pairs  per  millimeter. 

As  shown  in  Figure  3  the  dissector  system  with 
its  associated  control  electronics  is  interfaced  to 
the  PDP-S.E  minicomputer.  This  computer  buffers  a 
soar,  line  of  data  from  the  dissector  interiace  and 


FIGURE  3.  Image  Dissector  and  Digitization  (IDD) 
Facility 


writes  che  data  on  a  9-track,  computer  compatible 
tape  system,  one  line  at  a  time.  The  PDP-S/E  also 
has  the  ability  Co  perform  a  real-time  photonormal- 
ization  on  each  scan  line  of  data  for  removal  of  the 
nonuniform  response  characteristic  of  the  phoco¬ 
cathode  (typically  +11%  variation).  The  correction 
to  be  applied  is  automatically  determined  from  data 
collected  in  a  special  calibration  node  of  operation. 

The  images  are  digitized  to  8  bits,  giving  256 
image  intensity  levels.  An  alternate  circuit  is 
selectable  which  square  roots  the  video  prior  to 
digitization.  In  this  manner  256  levels  of  image 
amplitude  information  is  obtainable.  The  dynamic 
range  achievable  with  the  IDD  is  approximately  40 
dB,  compared  with  a  20  dB  dynamic  range  typically 
achievable  if  the  image  were  recorded  directly  on 
film. 

The  concrol  electronics  enable  the  user  to  se¬ 
lect  the  number  of  picture  elements  (pixels)  digi¬ 
tized  per  scan  line  up  to  a  maximum  of  3000  per  30 
ran  scan  length.  In  addition,  the  scan  length  can  be 
selected  by  setting  thumb  switches  to  the  desired 
scan  start  and  scan  stop  address  or  by  adjusting 
the  dissector  deflection  voltages.  The  integration 
time  per  pixel  is  selectable  in  increments  of  10 
Usees  from  10  usees  to  990  usees.  (Integration  is 
required  to  filter  out  electron  shot  noise  intro¬ 
duced  in  the  multiplier  stages.)  Finally,  a  point- 
skip  control  allows  the  user  to  vary  the  scan  line 
resolution  by  skipping  up  to  ten  pixels  between  di¬ 
gitizations.  Azimuth  sample  spacing  is  determined 
by  the  speed  of  the  precision  film  drive  used  to 
translate  the  input  film.  The  film  drive  unit  has 
a  long-term  stability  of  5  x  10  per  cent  and  a 
short-term  jitter  of  less  than  one  micron. 

The  IDD  system  is  capable  of  both  hardcopv  ar.c 
softcopy  displav  of  the  digitized  data.  Hardcopy  is 
provided  with  a  70  mm  CRT  film  recorder  while  soft- 
copy  is  accomplished  with  a  PEP-400  analog  scan 
converter  image  display  interfaced  to  a  PDP-11  40 
control  computer. 

The  hardcopy  CRT  film  recorder  is  interfaced  to 
the  PDP-8/E  and  used  for  displaying  digitized  data 
stored  on  the  magnetic  tape  unit.  The  CRT  is  a 
fiber  optics  direct  coupled  line  scan  recorder 
capable  of  recording  images  up  to  3000  elements  in 
range  and  virtually  any  length  in  azimuth  (limited 
by  film  length).  The  recording  spot  size  is  ap¬ 
proximately  20  microns  in  diameter  and  recording 
linearity  is  +0.1  per  cent  in  range  and  +  1  per 
cent  in  azimuth.  Although  not  as  accurate  as  the 
digitization  process,  this  display  technique  is 
useful  for  hardcopy  visual  presentation  of  digi¬ 
tized  imagery. 

The  softcopy  display  consists  of  a  PEP-^00 
Conrac  display  interfaced  to  a  PDP-11; -0  minicom¬ 
puter  which  receives  scan  dat3  via  the  PDP-S.  E  con¬ 
trol  computer.  The  PEP--.00  is  an  analog  scan  con¬ 
verter  which  utilizes  a  silicon  dioxide  image  sto¬ 
rage  tube.  The  storage  tube  is  capable  of  holding 
a  1024  x  1024  image  matrix  for  15  minute  before  it 
fades  and  has  available  32  shades  of  grey  for  dis¬ 
play.  The  Conrac  TV  monitor  is  a  1229  line  displav 
synchronized  to  the  PDP-400  controller.  In  addi¬ 
tion  to  the  conventional  contrast  ana  brightness 
control,  the  PEP-400  controller  has  a  zoom  control 
and  selective  X-Y  image  position  control  which  pro¬ 
vides  more  detailed  image  examination. 
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The  softcopy  display  provides  a  capability  for 
real-time  operator  interactive  processing.  Special 
interactive  processing  techniques  such  as  coherent 
optical  spatial  filtering  are  possible  on  a  real¬ 
time  basis.  The  principle  modes  of  operation  are: 

(1)  Direct  display  of  data  during  digitization; 

(2)  direct  display  of  Fourier  transformed  data 
during  digitization;  (3)  data  plots  of  FFTs  of 
image  data  for  initial  set  up  and  system  perfor¬ 
mance  evaluation;  and  (4)  data  tape  verification 
after  digitization. 

Examples  of  SAR  imagery  processed  and  digitized 
on  the  hybrid  IDD  facility  are  shown  in  Figure  4. 
The  images  were  produced  with  the  ERIM  four  channel 
multiplexed  radar  system  [5,6]  which  on  a  single 
pass  produces  phase  histories  at  two  wavelengths 
(3.2  and  23  cm)  with  two  polarization  configura¬ 
tions  for  each  wavelength  (horizontal  transmit, 
horizontal  receive  and  horizontal  transmit,  verti¬ 
cal  receive).  The  phase  histories  were  processed 


a)  3.2  cm  Horizontal-Horizontal  Polarization 


b)  3.2  era  Horizontal-Vertical  Polarization 


c)  23  cm  Horizontal-Horizontal  Polarization 


d)  23  cm  Horizontal-Vertical  Polarization 
FIGURE  -t.  Sar  Imagery  Optically  Processed,  Digi¬ 
tized  and  Displayed  using  the  Hybrid 
IDD  Facility. 


to  30  ft  resolution  in  both  range  and  azimuch  and 
were  digitized  with  6.67  ft  sample  spacing  in  boch 
coordinates.  The  images  were  then  reconstructed  on 
che  CRT  hardcopy  film  recorder.  The  area  imaged 
was  a  portion  of  the  NASA  soil  moisture  test  site 
near  Phoenix,  Arizona. 

DIGITAL  SAR  SIGNAL  PROCESSING  AND  IMAGE  EXPLOITA¬ 
TION 

Digitized  phase  histories  and  images,  including 
those  generated  on  the  hybrid  IDD  system  serve  as 
inputs  to  the  digital  radar  data  processing  and  ex¬ 
ploitation  facility.  The  goal  of  this  facility  is 
twofold:  (1)  the  development  of  two-dimensional 

digital  signal  processing  techniques  applied  to 
SAR  image  formation;  and  (2)  the  development  of 
manual  and  automatic  exploitation  techniques,  which 
will  improve  the  practicability  of  imaging  radars 
in  various  applications.  The  facility  comprises  a 
minicomputer  based  system  of  specialized  hardware 
in  a  multi-user  operating  environment.  The  system 
has  been  specifically  tailored  for  efficient  SAR 
processing  and  image  exploitation. 

A  block  diagram  of  the  digital  facility  appears 
in  Figure  5.  The  basic  component  of  the  system  is 
a  Digital  Equipment  Corporation  PDP-11/43  central 
processing  unit  which  utilizes  262,144  bytes  of 
on-line  memory  through  the  use  of  memory  manage¬ 
ment  hardware.  An  advanced  operating  system  (RSX- 
11D)  allows  multiple  real-time  tasks  to  run  con¬ 
currently  on  the  machine,  with  a  single  task  ad¬ 
dressing  up  to  65,536  bytes  of  memory.  Intertask 
communication  is  permitted,  which  implies  that 
several  tasks  could  potentially  be  working  on  a 
single  processing  problem.  Several  computer  ter¬ 
minals  are  provided  so  that  several  users  tar.  be 
utilizing  the  facility  simultaneously. 

Special  provisions  have  been  made  for  offline 
storage.  Due  to  the  vast  amounts  of  data  associ¬ 
ated  with  SAR  phase  histories  or  images,  a  116 
million  byte  random  access  disk  unit  has  been  in¬ 
cluded  in  addition  to  a  DEC  RK05  disk  unit  (2.4 
million  bytes)  used  for  operating  system  and  user 
program  storage.  The  magnetic  tape  system  consists 
of  two  9-track  and  one  7-track  dual  density  type 
drives,  allowing  tape-to-tape  processing  in  formats 
compatible  with  virtually  all  industry  standards, 
except  the  newer  6250  bpi  format. 

To  aid  in  implementing  SAR  digital  processing 
algorithms  a  hardwired  FFT  processor  has  been  in¬ 
cluded  in  the  svstem.  The  FFT  processor  (Time  Data 
Inc.  FPE4)  can  perform  a  1024  point  complex  FFT  in 
200  msec.  The  data  frame  size  for  the  FFT  is  con¬ 
trollable  from  4  complex  points  up  to  -096  complex 
points  in  powers  of  2.  In  addition  to  real  and 
complex,  direct  and  inverse  transforms,  the  proces¬ 
sor  can  perform  frequency  domain  Hanning  filtering 
as  well  as  auto  and  cross  spectrum  averaging.  Soft¬ 
ware  has  been  generated  which  uses  the  FFT  proeessei 
to  perform  a  large  array  two-dimensional  FFT  (up  to 
1024  x  1024  complex  points). 

.An  interactive  digital  image  displav  is  current¬ 
ly  being  integrated  into  the  digital  facility.  The 
RAMTEK  display,  which  has  a  solid-state  shift  re¬ 
gister  refresh  memory,  is  capable  of  displacing 
digital  imagery  in  two  modes  (1)  512  x  512  x  "  bits 
on  a  black  and  white  TV  monitor;  and  (2)  25n  x  25o 
on  a  color  TV  monitor  with  6  bits  each  controlling 
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che  red,  green  and  blue  intensities  for  each  pic¬ 
ture  element.  Overlay  channels  are  provided  so 
that  graphical  information  such  as  object  outlines 
or  annotation  data  can  be  superimposed  on  the  dis¬ 
played  images  without  destroying  the  contents  of 
the  display  memory.-  In  addition,  a  trackball  con¬ 
trolled  cursor  allows  the  user  to  extract  position¬ 
al  information  from  the  displayed  images.  A  dia¬ 
gram  of  the  interactive  display  system  is  shown  in 
Figure  6. 

An  additional  feature  of  the  display  are  the 
programmable  table-lookup  memories  installed  be¬ 
tween  the  display  memory  and  D/A  converters  which 
provide  the  video  intensities.  Under  program  con¬ 
trol  these  memories  can  be  loaded  with  intensity 
transfer  functions  which  are  then  applied  to  the 
images  on  a  point-by-point  basis  in  real  time  as 
they  are  displayed.  The  transfer  function  can  be 
constructed  to  provide  contrast  enhancement,  thres¬ 
holding,  image  intensity  inversion,  pseudo-color 
encoding  or  any  of  a  number  of  single  point  image 
manipulations,  the  results  of  which  are  nearly 
instantaneously  visible  Co  the  operator.  The  dis¬ 
play  will  be  used  to  view  digitally  processed  SAR 
imagery  and  to  examine  image  exploitation  techni¬ 
ques  . 
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FIGURE  5.  Digital  Radar  Data  Processing  and 
Exploitation  Facility 


SAR  digital  processing  consists  primarily  of 
two-dimensional  pulse  compression  techniques 
applied  to  digital  signal  histories.  The  digital 
facility  serves  as  a  general  purpose  cool  for 
examining  such  techniques.  The  2-D  pulse  compres¬ 
sion  is  generally  implemented  as  two  one-uire.ns  i  te¬ 
al  digital  matched  filtering  operations  in  con-unc¬ 
tion  with  range  and  azimuch  prefiltering  operations 
analogous  Co  the  frequency  plane  filtering  used  in 
the  optical  processing  technique  previously  de¬ 
scribed.  The  FFT  hardware  is  useful  for  performing 
che  various  types  of  digical  filtering  operations. 
The  digital  facility  has  been  structured  for  highly 
controlled  experimental  processing  raLher  than  em¬ 
phasizing  real-time  applications. 

Once  digital  SAR  imagery  is  available,  exploi¬ 
tation  processing  can  be  applied.  Such  techniques 
are  directed  towards  either  improving  manual  inter¬ 
pretation  or  implementing  machine  automated  inter¬ 
pretation  of  the  imagery.  The  digital  facility  is 
applicable  to  both  types  of  analysis. 

The  goal  of  processing  for  manual  interpreta¬ 
tion  is  to  increase  the  flow  of  information  from 
display  to  observer.  Most  conventional  image  pro¬ 
cessing  techniques  are  applicable  here,  ranging 
from  interactive  extraction  of  mensuration  data 
using  the  cursor,  to  spatial  frequency  filtering 
for  image  enhancement.  Emphasis  is  being  placed 
on  single  point  operations  to  take  full  advantage 
of  the  real-time,  programmable  enhancement  capa¬ 
bility  of  the  image  display. 

Automated  interpretation  involves  the  use  of 
pattern  recognition  algorithms  on  features  ex¬ 
tracted  from  the  radar  imagery.  Applications  in¬ 
clude  automatic  target  identification  and  classi¬ 
fication  of  land  use  or  conditions  on  the  basis  of 
multi-channel  radar  returns.  For  examt  le  ,  the  : 
images  shown  in  Figure  4  are  being  analyzed  in  an 
effort  to  determine  soil  moisture  conce-r ,  surface 
roughness  characteristics,  and  vegetation  type 
associated  with  the  field  patterns  in  the  image. 

The  four  channels  have  been  digitally  registered 
and  multivariate  discriminate  analysis  is  being 
applied  to  derive  decision  rules  based  on  tn.e  f  -ur 
radar  returns  for  each  picture  element  (*].  In 
addition,  texture  analysis  [ 8 ]  is  being  applied  to 
individual  images  in  an  effort  to  determine  the 
relationship  between  each  surface  condition  and 
the  spatial  characterist ics  of  the  radar  returns. 


SUMMARY 

A  system  has  been  described  which  is  desi  one-.: 
specifically  tor  signal  processing  and  image  ex¬ 
ploitation  of  SAR  data.  The  total  facility  uti¬ 
lizes  coherent  optical,  hybrid  opt ica  1- a  i  g  i : a  I  . 
and  wholly  digital  processing  technologies  to 
generate  and  effectively  display  SAR  imagery.  The 
optical  system  uses  a  state-o  f-the-  irt  tilted  xAo, 
processor  to  process  optically  recorded  S.AK  signal 
histories  into  imagery.  The  hybrid  facility  maxes 
use  of  an  image  dissector  device  as  a  bridge  be¬ 
tween  the  optical-analog  and  digital  realms  to 
effectively  combine  the  speed  and  power  of  coherent 
optical  processing  with  the  flexibility  and  rro  :- 
sion  ot  digital  cechniqties.  A  unique  digital  s-  s- 
ten.  which  includes  a  hardwired  FFT  process  -  in.-, 

art  interactive  digit ai  image  display  can  -non 
used  to  investigate  digital  SAR  processing  ox:  -.  x, 
manual  and  automated  interpret.it  ion  terhni:  ■. 
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FIGURE  6.  RAMTEK  Digital  Image  Display  Block  Diagram 
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APPENDIX  D 
DATA  TAPE  FORMAT 

The  following  is  a  description  of  ERIM  digital  data  tape  format 
(SAR  IPL  data): 

1.  Each  CCT  (Computer  Compatible  Tape)  has  been  recorded  on  9 
track  tape  at  a  density  of  6250  6PI. 

2.  There  are  8  files  on  each  tape,  with  2  end  of  volume  marks 
immediately  following  the  last  file. 

3.  The  first  record  of  each  file  contains  a  header  record  from 
the  original  data.  The  first  96  bytes  of  this  record  are 
used  to  store  the  following  information: 

A.  Data  type  number 

B.  File  number 

C.  Flight  date 

D.  Band 

E.  Polarization 

F.  Bite 

G.  File  size  (prior  to  geometric  correction) 

H.  Resolution  in  meters  (range) 

I.  Data  Compression  method  (square  root) 

J.  Processing  date 

All  header  data  is  in  an  ASCII  Convention. 

4.  Processed  SAR  data  (beginning  with  record  2  of  each  file) 
consist  of  fixed  length  unformatted  records.  The  SAR  data 
is  in  a  byte  convention  (quantized  to  3  bits,  one  byte  per 
image  element).  Information  on  the  number  of  records  and 
record  length  in  bytes  of  each  file  has  been  included,  as 
well  as  an  example  of  header  information. 
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